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A  detailed  knowledge  of  the  gap  structure  for  the  Fe-pnictide  superconductors  is  still  rather  rudimentary,  with 
several  conflicting  reports  of  either  nodes,  deep  gap  minima,  or  fully  isotropic  gaps  on  the  Fermi  surface  sheets, 
both  in  the  kx  —  ky  plane  and  along  the  c— axis.  In  this  paper  we  present  considerations  for  electronic  Raman 
scattering  which  can  help  clarify  the  gap  structure  and  topology  using  different  light  scattering  geometries. 
Using  density  functional  calculations  for  the  Raman  vertices,  it  is  shown  that  the  location  of  the  gap  minima 
may  occur  on  loops  stretching  over  a  portion  of  the  c-axis  in  BafFeq  94(300.06)2- 


PACS  numbers:  74.25.nd,74.70.Xa,74.20.Pq,71.15.Mb 

Since  the  discovery  of  high  temperature  superconductivity 
in  the  iron  pnictides,  identifying  the  structure  of  the  super¬ 
conducting  order  parameter  has  remained  a  tantalizing  and 
elusive  problem.  Knowledge  of  the  momentum  dependence 
is  a  major  step  towards  identifying  the  interactions  that  drive 
the  superconductivity.  A  variety  of  experimental  techniques 
probing  the  low  energy  properties  of  the  pnictides  in  the  su¬ 
perconducting  state  have  given  contradictory  information  on 
the  angular  dependence  of  the  superconducting  gap  around 
the  Fermi  surfaces  (FS).  Angle-resolved  photoemission  on 
NdFeAsOo.9Fo.i 1,  K-doped2  and  Co-doped3  BaFe2As2  have 
identified  rather  isotropic  superconducting  gaps  on  both  the 
hole  and  electron  pockets  in  the  Brillouin  zone  (BZ).  On 
the  other  hand,  transport  and  thermodynamic  measurements, 
which  do  not  have  angular  resolution,  have  indicated  the  pres¬ 
ence  of  nodes  in  some  materials  and  small  but  finite  gaps 
in  others4.  Recently,  angle-resolved  thermal  conductivity  in 
BaFe2(As,P)2  have  indicated  large  gap  variations  around  the 
Fermi  surface,  with  gap  nodes  or  deep  minima  located  at  45° 
with  respect  to  the  principal  BZ  axes5. 

Electronic  Raman  scattering  experiments  indicated  very 
small  gaps  located  on  the  electron  Fermi  surfaces  in 
Ba(Feo.94Coo.o6)26-  However,  due  to  the  relative  insensitiv¬ 
ity  of  the  Raman  vertices  to  momentum  in  the  canonical  tight 
binding  (TB)  scheme  the  location  of  the  nodes  could  not  be 
determined  uniquely.  The  situations  is  additionally  obscured 
by  the  fact  that  the  electronic  structure  of  the  122  materials 
is  inherently  more  3D  than  that  of  the  1111  pnictides,  and 
simplified  2D  tight  binding  models  can  be  rather  misleading7. 
Moreover,  the  crystallographic  symmetry  of  the  122  structure 
is  body-centered  tetragonal  (bet),  M/mmm ,  therefore  folding 
down  the  Fermi  surfaces  from  the  one-Fe  Brillouin  zone  to  the 
two-Fe  zone  is  less  trivial  than  in  the  1111  structure.  Namely, 
while  in  the  1111  structure  the  electronic  pockets  at  the  X 
point  and  Y  point  in  the  unfolded  zone  are  related  by  a  simple 
7f/2  rotation,  so  that  after  folding  they  form  e-pockets  with  the 
full  4-fold  symmetry,  in  the  122  structure  the  symmetry  oper¬ 
ation  connecting  the  two  unfolded  pockets  is  the  same  rotation 
plus  a  shift  along  kz  by  n/c.  Thus  for  all  kz  but  kz  =  n/2c  the 
folded  FS  cuts  have  only  a  2-fold  symmetry. 


Incidentally,  for  the  same  reason  the  standard  method  of 
constructing  a  TB  Hamiltonian  by  transforming  the  Bloch 
functions  into  the  Wannier  functions  leads  to  a  TB  model  that 
attempts  to  describe  the  outer  barrel  as  a  single  band  in  the 
unfolded  BZ7,  while  in  reality  it  is  formed  by  anticrossing  the 
X  and  the  Y  pockets.  In  other  words,  such  TB  models  feature 
unphysically  and  artificially  large  interlayer  hoppings  (some 
interlayer  hoppings  in  Ref.  7  appear  to  be  more  than  50% 
of  the  nearest-neighbor  in-plane  hopping),  while  in  reality  the 
scale  of  the  kz  dispersion  is  set  by  the  hybridization  gaps,  and 
by  the  difference  in  size  (but  not  orientation)  between  the  up¬ 
per  and  the  lower  panel  in  the  Fig.  4  of  Ref.  7  . 

For  this  reason,  in  the  122  systems  extracting  quantitative 
information  from  the  experiment  using  an  infolded  2D  tight 
binding  model  becomes  questionable,  and  working  in  a  folded 
3D  TB  model  is  not  much  simpler  than  straight  off  first  prin¬ 
ciple  calculations.  Therefore,  in  order  to  gain  more  accurate 
information  about  the  location  of  gap  minima  on  the  Fermi 
surface  we  have  performed  all-electron  band  structure  calcu¬ 
lations,  as  described  in  Ref.  8  and  numerically  computed  the 
relevant  Raman  vertices.  We  have  confirmed  the  result  ob¬ 
tained  qualitatively  by  Muschler  et  al.6  that  two  of  the  three 
principal  Raman  symmetries,  A\g  and  lh:z  (all  notations  in  the 
crystallographic  unit  cell),  probe,  respectively,  an  average  gap 
and  the  gap  on  the  electron  Fermi  surface,  and  that  for  the  third 
symmetry,  B  \  „ .  the  matrix  elements  are  mostly  small  on  all  the 
FSs.  However,  we  found  that  on  the  outer  electronic  FS  there 
exist  “hot  spots”  where  this  matrix  element  is  sizeable,  a  fact 
that  could  not  be  seen  in  unfolded  TB  calculations.  Therefore 
we  can  not  only  distinguish  between  the  hole  and  the  electron 
FSs  in  terms  of  the  gap  minima,  but  can  also  place  them  at  or 
near  these  hot  spots  (which  appear  near  the  [110]  direction). 

First,  we  discuss  the  so-called  “effective  mass  approxima¬ 
tion”  as  applied  to  this  compound.  The  inelastic  light  (Raman) 
scattering  cross  section  in  the  limit  of  small  momentum  trans¬ 
fers  is  given  in  terms  of  correlation  functions  of  an  effective 
charge  density9,  p  =  L„,k  a  7«(k,  (Oi:  filj)c(:!  k  (Tc„^k.cT, where 

cnko  O'n.k.rr)  creates  (removes)  an  electron  with  momentum  k 
and  spin  cr  from  band  n.  In  principle,  y  depends  both  on  light 
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FIG.  1 :  (Color  online )  Raman  response  for  xx  and  xy  polarization  for 
three  different  excitation  lines  as  indicated.  All  spectra  have  the  Bose 
factor  removed.  Apart  from  correcting  for  the  instrumental  response, 
the  overall  intensity  is  as  measured  and  has  not  been  shifted. 


scattering  polarization  geometries  e'^  and  the  frequencies  of 
incoming  (scattered)  light  (Ol(s\  In  the  limit  where  ft) 'W  are 
much  smaller  than  any  relevant  interband  transition  frequency, 
the  scattering  amplitude  simplifies  into  the  well-known  effec¬ 
tive  mass  tensor 


7n(k,  (Oi,s  — >■  0)  = 


L  * 

a,[S=x,y,z 


i  d2e„(k ) 
adkadkp 


(1) 


For  tetragonal  materials,  three  different  scattering  geometries 
can  be  classified  according  to  the  point  group  transformation 
properties.  These  are  given  as 

Aig  :  Jn (k)  =  [d2£n(k)/dk2  +  d2£„(k)/dk2]/2 
Big  :  7„(k)  =  [d2en(k)/dk;-d2en(k)/dk2\/2 
B2g :  7«(k)  =  d2en(k)/dkxdky.  (2) 


Thus  the  curvature  of  the  bands  and  the  light  scattering  po¬ 
larization  orientations  determine  which  carriers  are  involved 
in  light  scattering  in  different  bands  and  regions  of  the  BZ. 
While  not  having  the  definitive  angular  resolution  as  photoe¬ 
mission  to  probe  single  particle  properties,  Raman  scattering 
can  selectively  project  effective  charge  density  fluctuations  in 
different  regions  of  the  BZ,  uniquely  providing  k— projected 
dynamical  two-particle  properties. 

In  most  systems,  such  as  the  cuprates  and  pnictides,  inci¬ 
dent  light  energies  are  typically  high  enough  to  directly  ex¬ 
cite  interband  charge  transfer  and  thus  the  effective  mass  ap¬ 
proximation  may  be  not  valid  for  these  systems.  When  this 
happens,  strong  resonance  effects  occur  and  the  intensity  of 
the  electronic  Raman  scattering  varies  greatly  with  the  inci¬ 
dent  photon  energy,  as  for  instance  in  cuprates.10  However, 
as  shown  in  Fig.  1  in  the  pnictides  the  resonance  effects  are 
weak  and  thus  off-resonance  conditions  for  the  scattering  am¬ 
plitude  applies.  Therefore  we  assume  that  the  effective  mass 
approach  holds  for  the  pnictides  and  examine  the  full  band 
and  momentum  dependence  of  the  scattering  amplitudes  as  a 
way  of  pinpointing  the  gap  extrema  structure. 


In  the  superconducting  state,  the  Raman  response  for  B\g 
and  B2g  symmetries  is  given  as  a  weighted  average  of  the 
Tsuneto-Maki  function  A„  over  the  FS  for  each  band  nn, 
Xy,y(Q)  =  .k  7?(k)A„(k,fI),  with  the  Raman  vertices  y 
given  in  Eq.  2.  The  Tsuneto  function  is 


A,,(k,n)  =  tanh 


E„(k)  4  |  A„(k)  |2  /En(k) 

2 kBT  4E2(k)-(hQ.  +  i8 )2’ 


(3) 


where  A„(k)  is  the  superconducting  gap  in  band  n ,  having 
quasiparticle  energy  dispersion  £„(k),  and  E2(k)  =  £2(k)  + 
A“(k).  Enforce  particle  number  conservation,  one  gets 


X(&)  =  *r,y(Q)  -  Xy,  l(&)Xl,y(&)/Xl,l(&),  (4) 

with  =Xuy,(^)  =E«,k7n(k)^(k>^)  and  Xi,i  (ki)  = 

En,k^n(k ,kl).  Therefore,  while  for  both  B\g  and  B2g  each 
contribution  to  the  Raman  vertex  can  be  separated  into  con¬ 
tributions  from  single  bands,  the  A  i?  contributions  mix  mass 
fluctuations  in  each  band  and  between  bands. 

In  Ref.  6,  the  interplay  of  the  energy  gap  and  the  symmetry 
dependence  of  the  Raman  vertices  was  considered,  where  it 
was  argued  that  B2g  polarizations  largely  sampled  the  energy 
gap  on  the  electron  bands,  while  A  was  more  sensitive  to  the 
hole  bands  in  the  pnictides.  Here  we  explore  the  form  for  the 
Raman  vertices  in  more  detail  using  LDA. 

The  Raman  vertices  shown  in  Eq.  (2)  are  directly  com¬ 
putable  from  the  one-electron  band  structure.  We  have  per¬ 
formed  full-potential  LAPW  calculations  using  the  WIEN2k 
package  and  an  experimental  crystal  structure,  including  the 
spin-orbit  interaction.  The  doping  is  accounted  for  in  the 
virtual  crystal  approximation.  The  derivatives  of  the  eigen¬ 
states  were  calculated  numerically  using  a  31x31x31  /.'-point 
mesh  in  the  bet  BZ.  We  neglect  any  momentum  interplay  be¬ 
tween  the  Raman  vertices  and  the  superconducting  gaps  in 
the  Tsuneto-Maki  function,  and  examine  just  the  character 
of  the  Raman  vertices  on  each  band  for  each  polarization 
geometry. 12  The  calculated  quantities  are  derived  from  eval¬ 
uating  Eq.  4  as 


*£y  =  Itf(k)5Mk)] 

n:  k 


(L,,k7n(k)g[en(k)])2 

L,.k<S[£„(k)] 


(5) 


This  can  be  simplified  by  writing  y„(k)  =  Ay,(k)  +  y„,  with 
%  =  Lk7«(k)^[£«(k)]/Lk^[e»(k)]  the  average  of  the  Raman 
vertex  for  band  n  over  Fermi  surface  sheet  n.  The  response 
can  be  separated  into  an  intraband  piece,  which  is  the  only 
contribution  for  B\g  and  B2g, 

Xintraband  =  Ay,  (  k )  ~  8  \£fl  ( k )] ,  (6) 

n.k 


and  an  interband  contribution  only  appearing  for  A  |  „  channels 
'LnNnik  ~  {Y.nNn7n)2 /Y.nNn-  The  last  part  can  be  written  as  a 
band-dependent  sum: 

Xinterband  =  irrYn  T" '  NmiYu  Ym )  Xu -ini (‘rband ■  (2) 
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with  Nt  =  Nn  the  total  DOS.  The  results  are  shown  in  Table 
I  and  are  defined  as  follows: 

M,  =  £S[e„(k)] 

k 

Yn.ap  =  ^[<32e«(k)/<9^a^]5[e„(k)]. 

k 

Y„,  apyS  =  E  [^2e«  (tt/dkadkp][d2en  (k) /d  k7dks\ 8  [en  (k)] , 

k 

such  that  the  B\g  polarizations  select  (jn.xxxx  —  Yn,xxyy) /2,  Big 
selects  Yn,xyxy  The  unscreened  intraband  part  for  A\g  is  given 
by  (Yn, vax  +  7;i,x*yv)/2-  which  is  reduced  with  screening  to 

(Yn,xxxx  T  Yn.xxyy) /2  Nn y~ ■ 


hi 

h2 

h3 

el 

e2 

Nn 

10.85 

16.4 

12.97 

12.33 

10.84 

( Yxx  +  Yyy)/^ 

-18.2 

-9.9 

-12.3 

9.1 

19.7 

Y 

-1.7 

-0.6 

-1.0 

0.74 

1.8 

(Yxxxx  +  Yxxyy)/2 

39.2 

9.4 

15.9 

19.7 

47.8 

A\g  (intraband) 

7.8 

2.6 

2.9 

15.1 

12.7 

A\g  (interband) 

28.6 

3.6 

7.7 

8.6 

38.9 

Big 

2.2 

4.3 

6.3 

32.8 

8.5 

B2g  (averaged) 

4.8 

8.5 

5.2 

21.8 

13.8 

TABLE  1:  Values  of  the  Raman  vertices  and  related  quantities  in¬ 
cluding  spin-orbit  coupling  around  the  hole  (hl,h2,h3)  and  electron 
(el,e2)  Fermi  surfaces,  on  a  per  cell  basis,  in  Ry  units  (energy  in 
Ry,  distance  in  Bohr  radii,  momentum  in  inverse  Bohr  radii).  Here 
hi. ..  h3/el.  ..e2  denote  the  outermost.. .innermost  hole/electron  barrel 
Fermi  surface,  centered  at  TIM  in  the  Brillouin  zone,  as  indicated  in 
Fig.  2. 

The  LDA  results  shown  in  Table  I  indicate  that  both  B\g 
and  Z?2g  polarizations  give  vertices  which  are  dominated  by 
contributions  from  the  electron  FS,  in  agreement  with  the  as¬ 
sessment  of  Ref.  6.  Moreover  we  note  that  B\g  vertices  are 
specifically  large  only  on  a  single  electron  FS  sheet,  given 
by  the  outer  FS  pocket  centered  at  M.  For  the  A  |  g  contri¬ 
bution,  while  the  unscreened  Raman  vertex  Yxxxx  +  Yxxvy  has 
large  contributions  from  the  hole  band  hi  and  electron  band 
e2,  backflow  largely  cancels  the  intraband  contribution,  leav¬ 
ing  a  substantial  interband  contribution  deriving  largely  from 
these  same  bands.  These  results  support  the  conjecture  based 
on  symmetry  in  Ref.  6  that  the  hole  bands  would  be  probed 
for  A  ij,  configurations  more  so  than  for  B\g  and  Bjy- 

The  LDA  results  indicate  that  the  A{g  interband  portion  in¬ 
volves  also  a  strong  contribution  from  scattering  between  the 
electron  band  e2  and  the  hole  band  hi,  coming  from  different 
signs  of  the  Raman  vertices  in  the  two  bands.  As  discussed 
in  Ref.  1 1  the  divergence  at  twice  the  gap  edge  which  may 
appear  in  crossed  polarization  orientations  (such  as  Big.2g)  is 
removed  in  single  band  superconductors  for  A\g  symmetries 
due  to  the  gauge  invariant  backflow  which  screens  uniform 
charge  displacements  over  many  unit  cells.  In  multiband  sys¬ 
tems,  such  as  the  pnictides,  the  superconducting  response  for 


FIG.  2:  (Color  online)  Left:  Raman  vertices  for  B\g  symmetry,  eval¬ 
uated  around  each  Fermi  surface  sheet  projected  onto  the  x  —  y  plane, 
that  is,  looking  down  the  k,-direction.  The  parts  of  the  Fermi  sur¬ 
face  where  the  value  of  B\g  Raman  vertex  is  larger  than  the  average 
are  shaded  red.  The  gold  ellipses  mark  the  area  where  the  deep  gap 
minima  or  nodes  are  expected.  Right  top:  The  same,  viewed  along 
the  M  —  T  line,  as  shown  by  the  arrow  in  the  left  panel;  only  the 
outer  electron  FS  pockets  are  shown.  Right  bottom:  The  same  view, 
but  now  the  parts  where  the  Fermi  velocity  is  larger  than  60%  of  its 
maximal  value  are  red. 

A \g  symmetries  may  be  written  as11 

(fl)  =  I  (Yn  (k)  -  ?m2  K  (k,  fl) ,  (8) 

72,  k 

with  A,,  given  by  Eq.  (3).  The  Tsuneto-Maki  function  A„  pos¬ 
sesses  divergences  when  fi  =  2  |A„(k)|.  y(D)  is  another  way 
of  representing  the  frequency-dependent  backflow  in  the  su¬ 
perconducting  state  for  A\g  symmetries: 

?(fl)  =  L7»(k)^(k)fl)/£AB(k)fl).  (9) 

72, k  77, k 

Each  separate  divergence  occuring  at  the  gap  edge  in  each 
band  in  A,,(k.£2)  weighs  exactly  those  momentum  points  in 
y„  (k)  for  which  the  divergence  occurs,  and  via  Eq.  8  the  di¬ 
vergence  is  cancelled.  This  holds  for  each  divergence  coming 
from  each  band  n,  and  therefore  we  conclude  that  no  diver¬ 
gence  may  appear  in  A\g  symmetries  for  multiband  systems. 
This  conclusion  holds  only  for  Cooper  pair  creation  via  light 
scattering  but  does  not  consider  additional  divergences  com¬ 
ing  from  charge  or  spin  collective  modes.13,14 

This  predicted  absence  of  an  A\g  divergence  is  in  general 
agreement  in  the  observed  data  on  the  pnictides,  where  no 
clear  2A  scattering  peak  is  seen  for  A  ig  scattering  geome¬ 
tries.  While  sharp  collective  modes  may  appear  due  to  the 
oscillation  of  pairing  amplitudes  in  the  different  bands  for  A\g 
symmetries15,  no  clear  indication  of  such  a  mode  has  been 
found.  Similarly,  although  collective  modes  have  been  pre¬ 
dicted  to  also  occur  in  B2g  symmetries  due  to  final  state  exci- 
tons  having  orthogonal  symmetry  to  the  superconducting  pair 
state16,  it  is  still  an  open  question  whether  these  modes  can  be 
uncovered  in  the  data. 

The  fact  that  in  the  experiment  hardly  any  effect  of  super¬ 
conductivity  is  observed  in  the  B\g  polarization,  while  the  cal¬ 
culations  predict  a  noticeable  coupling  with  the  outer  elec¬ 
tronic  FS,  suggests  that  this  coupling  is  not  uniform,  but  is 
small  wherever  the  gap  is  sizeable,  and  where  the  coupling 


4 


is  large  the  gap  is  small.  In  order  to  test  this  hypothesis,  we 
have  plotted  the  effective  mass  asymmetry,  appearing  in  B\g 
polarizations  as  |<?2£„(k )/dkx  —  <?2£„(k)/clk2|,  in  Fig.  2.  We 
observe  that,  first,  this  difference  is  small  on  all  FSs  except  the 
outer  e-barrel,  and,  second,  it  is  also  small  on  the  most  of  that 
FS  as  well,  except  on  a  small  spot  located  at  an  angle  varying 
from  10  to  15°,  depending  on  kz.  As  Fig.  2  shows,  this  is  con¬ 
sistent  with  gap  nodes  or  deep  gap  minima,  forming  elliptical 
loops  on  the  outer  electron  FS  centered  around  the  | ,  0) 
points,  where  the  T  —  M  line  intersects  the  electron  pockets 
(indicated  by  the  ellipse  in  Fig.  2).  This  can  be  compared  with 
model  calculations  by  Graser  et  al ?  who  found  that  the  pair¬ 
ing  interaction  is  strong  in  the  Z  plane  (kz  =  n/c)  plane  and 
weaker  in  the  F  plane.  In  their  calculations  this  led  to  a  nodal 
loops  similar  to  ours,  but  located  at  the  hole  FS.  One  should 
be,  however,  careful,  when  making  a  quantitative,  as  opposed 
to  qualitative,  comparison,  with  these  calculations,  because, 
as  discussed  above,  Graser  et  al.  worked  in  an  unfolded  BZ 
wher  the  folding-induced  3D  dispersion  was  emulated  through 
enhance  kz  hopping  parameters. 

It  is  also  consistent  with  the  thermal  conductivity  data5 
(collected  on  another  122  material),  which  indicate  that  there 
are  gap  nodes  present  in  the  system,  yet  the  specific  heat  does 
not  seem  to  be  affected  by  them.  These  data  have  been  in¬ 
terpreted  in  terms  of  ultraheavy  hole  bands  that  contribute  at 
least  80%  to  the  DOS,  and  ultralight  electron  bands  that  dom¬ 
inate  the  transport.  However,  the  former  assumption  seems 
a  bit  too  strong  to  be  realistic.  Yet  if  in  BaFe2(As,P)2  the 
gap  structure  resembles  the  one  that  we  propose  here  for 
BaFei  88C00.12AS2,  so  that  the  gap  nodes  exist  only  for  a  finite 
range  of  kz.  the  rest  of  the  electronic  bands,  fully  gapped,  add 
their  contribution  to  the  hole  band  and  together  provide  80% 
or  more  of  the  total  specific  heat.  In  fact,  incidentally,  the 
Fermi  velocity  is  the  largest  close  to  the  “hot  spots”  that  we 
have  identified  as  locations  of  gap  nodes,  in  agreement  with 
the  idea  that  the  transport  is  dominated  by  near-nodal  parts  of 
the  FS  (Fig.  2). 

One  can  ask  a  question,  why  most  of  the  Fermi  surface, 
except  for  relatively  small  “hot  spots”,  shows  such  weak  cou¬ 
pling  in  the  B\g  channel.  This  is  actually  easy  to  address.  Let 
us  start  with  a  simple  model,  a  2D  parabolic  electronic  band. 
By  symmetry,  the  principal  axes  of  the  corresponding  ellipti¬ 
cal  FS  are  along  (11)  and  (11).  For  the  moment,  we  will  use 
the  center  of  this  FS  as  the  origin  and  the  principal  axes  as 
coordinates.  Then  £12  =  +  H2.\ k2.  The  Raman  vertex 

in  question  appears  like  a  ZLg  vertex  in  this  basis  and  is  pro¬ 
portional  to  d2en(kx,ky)/dkxdky  =  0.  Thus  in  this  model  this 
band  does  not  couple  in  the  B\g  channel  at  all. 

However,  this  is  not  necessarily  true  if  one  takes  into  ac¬ 
count  the  downfolding  of  the  electron  pockets  in  the  crystal¬ 
lographic  unit  cell.  After  the  downfolding  the  two  FSs  simply 
intersect.  In  our  model,  as  well  as  in  the  real  1 1  and  1111 
crystallography,  they  cross  along  the  kx  -  0  and  ky  =  0  di¬ 
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rections,  and,  because  they  have  different  parities,  there  is  no 
hybridization  between  them  and  the  newly  formed  degener¬ 
acy  is  not  lifted  (although  it  would  be  lifted  by  the  spin-orbit 
coupling).  On  the  other  hand,  in  the  122  symmetry  the  inter¬ 
sections  occur  at  general  points  in  the  Brillouin  zone  so  that  no 
symmetry  prevents  them  from  hybridization,  thus  forming  an 
inner  barrel  and  an  outer  barrel  (el  and  e2)  even  without  the 
spin-orbit.  To  gain  more  insight  into  this  problem  we  intro¬ 
duce  hybridization  into  our  model,  and  write  the  Hamiltonian 
near  the  anticrossing  point  as  Hu  =  t'  \  (k),H22  =  £2(k)-//i2  = 
H21  =  V,  where  V  is  a  hybridization  matrix  element.  The  in¬ 
verse  masses  of  the  resulting  two  bands  are 

d2S±  8V2|(/Ji  —  ll2)1kxky\ 

dkxdky  [(Ml-M2)2fe2-k2)2  +  4V2]3/2- 

This  is  essentially  an  approximate  expression  of  a  delta  func¬ 
tion  centered  at  kx  ±  ky  =  0  with  the  width  V  and  strength 
|/li  —  JU.2 1  ■  Thus  in  this  simple  model  we  have  vanishing  cou¬ 
pling  everywhere  except  hot  spots  where  the  bands  would 
cross  if  the  As  sublattice  would  not  have  a  lower  symmetry 
than  the  Fe  one. 

This  simplified  picture  only  qualitatively  applies  to  Bal22. 
The  bands  there  are  not  exactly  parabolic,  and,  most  impor¬ 
tantly  they  cross  not  at  the  kx,  ky  =  0  lines,  but  at  general  points 
that  can  be  considerably  removed  from  these  line,  depending 
on  kz.  Yet  our  model  demonstrates  qualitatively  how  the  “hot 
spots"  for  this  polarization  appear  on  the  FS,  and  why  they  are 
so  small  in  size. 

In  summary,  in  this  paper  we  have  paired  experimental  data 
on  the  electronic  Raman  scattering  in  the  superconducting 
state  of  Ba(Feo.94Coo.o6)2  with  the  first  principle  calculations 
of  the  effective  mass  fluctuations.  We  found  that  a  2D  TB 
model  in  the  unfolded  Brillouin  zone  is  not  sufficient  to  cap¬ 
ture  all  the  mass  fluctuation.  The  reason  is  that  because  of  the 
body-centered  symmetry  of  the  corresponding  crystal  struc¬ 
ture  downfolding  the  electron  bands  results  in  hybridizations 
and  substantial  band  anticrossings,  which  in  turns  produces 
large  mass  fluctuation  in  a  small  part  of  the  Fermi  surface. 
Our  results  suggest  that  closed  nodal  loops  (rather  than  verti¬ 
cal  or  horizontal  nodal  lines)  exist  in  this  compound  near  the 
intersection  of  the  outer  electron  Fermi  surface  with  the  T— M 
line.  It  is  likely  that  the  size  and  the  very  existence  of  these 
loops  is  doping-dependent,  but  this  is  beyond  the  scope  of  this 
paper. 
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